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ABSTRACT 

The Galactic bulge shows a double red clump in sight lines with |6| > 5° and 
—3° ^ ^ ^ 4°. It is interpreted as a signature of an X-shaped structure seen almost 
edge-on. We measure the proper motions of the stars belonging to the closer and 
the further arm of the X-shaped structure. The intrinsic kinematic properties of 
the two arms are found by incorporating information taken from the luminosity 
function. At b = —5° we find that the proper motion difference between the two 
arms is a linear function of the Galactic longitude for —0.1° < / < 0.5°, which 
we interpret as a streaming motion of the stars within the X-shaped structure. 
Such a streaming motion has not previously been reported. The proper motion 
difference is constant for —0.8° < / < —0.1°, which gives us an estimate of bulge 
rotation speed of 87.9 ± 8.2 km s^^ kpc^^. 

Subject headings: Galaxy: structure — Galaxy: kinematics and dynamics — stars: 
kinematics — astrometry 
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Introduction 



The bulge of the Galaxy is the closest such structure and the only one in which stars 
can be well resolved and studied individually. This in turn allows detailed investigation 
of properties of stellar populations within the bulge. T he Galactic bulge was not formed 



via ma 



or mergers 



found ( JNataf et al. 



3ut was developed via bar bucklin g (IShen et al. 



2010 : 



McWilham fc Zoccali 



20101). It was recently 



20101 ) that the Galactic bulge contains an 



X-shaped structure. This structure is aligned with a bar, which makes ~ 30° angle with 
Sun - Galactic center direction, and its arms are significantly extended along the line of 
sight direction. The color-magnitude diagrams (CMDs) constructed in some bulge areas 
show a double-peaked red clump (RC) because the line of sight crosses two arms of the 
X-shaped structure. In fields close to the Galactic disk the two RCs merge and appear 
single-peaked. These include the well studied low-extinction Baade window at Galactic 
coordinates {l,b) = (1°, —3.9°). The RC is also single-peaked in the areas where only one 



arm is seen. X-shaped structures are known in other galaxies (jWhitmore fc Bell 



19881) and 



may be centered (four arms cross at the center of the galaxy) or off-cei itered (arms cross i n 



two areas within the disk that do not coincide with the Galactic center; 



Bureau et al. 



20061 ). 



The double RC in the Galact i c bulg e was confirmed 



infrared data by ISaito et al. 



(12011 



20121 ). More recently 



using 2MASS an d VVV near 



Cao et al. 



(120131 ) analyzed the 



spatial density of the RC stars and found a clear signature of the X-shaped structure. 



Li fc ShenI ( 12012 



Shen et al. 



found properties of the X-shaped structure in the bulge model of 



( 120101 ) to be consistent with observations. 



Until now there have been relatively few studies devoted to the kinematics of both 
RCs. In some cases they were conducted at fields where the X-shaped structure is present 
but hard to be distinguished i.e., —5° ^ b < —3°. Three such e xamples wer e repo rted by 



Mc William fc Zoccali (i2010i) . The catalog of proper motions by iSumi et al. 



( I2004h based 
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on the second phase of the Optical Gravitational Lensing Experiment (OGLE-II) data 
was used to calculate a proper motion difference of 1.0 ± 0.06 milliarcseconds per year 
(mas yr~^) between the two arms of the X-shaped structure at the Baade window. The 



Vieira et al. 



catalo g of proper motions in the Plant field (/, 6) = (0°, —8°) compiled by 
( 120071 ) resulted in a proper motion d ifference in longitude o f 0.5 1 ± 0.18 mas yr~^, and 
of 0.19 ± 0.19 mas yr~^ in latitude. iMc William fc Zoccalij ( i2010i ) claimed also a SAa 
difference in the latitudinal proper motion dispersion without providing the figures 
themselves. These authors also recalled the radial velocity (RV) dif ferences between the 



faint and bright part of the RC obtained by 



Rangwala et al. 



(120091): -40 ± 11 km s'^ 



at (/, 6) = (5.5°, —3.5°), —4 ± 15 km s ^ a t the Baade window, a nd — 32 ± 11 km s ^ at 



(/,6) 



-5.0°, -3.5°). As pointed out by 



McWiUiam fc Zoccalil (12010[ l the results in the 



two fields with / p^ ±5° have magnitudes that are very close to the estimates of the model 

I fl [ 

predictions by iMao fc Paczvhskil f|2002i) bu t have the opposite sign. The Baade window 



was also studied by 



Babusiaux et al. 



(l2010f ) . An RV difference of 70 ± 30 km s ^ was founc 



after e xcluding metal-poor stars. The pro per motion differen ces found by 



Babusiaux et al. 



( I2OIOI ) in the Baade window based on the 



Sumi et al 



(|2004[ ) catalog were not statistically 



significant. 



De Propris et al. 



(120111 ) measured a nonsignifica n t RV difference of 12 ± 10 km s ^ 



between two RCs at the Plant window. 



Ness et al 



(120121 ) measured the RV offset 



of -30 ± 12 km s"^ at (/, 6) = (0°, -5°) and 7 ± 9 km s"^ in joint data for two 
fields: {l,b) = (0°, -7.5°) and (lb) = (0° -10°). Both figures were consistent with 



model predictions. lUttenthaler et al. 



(I2OI2I) found 4.4 ± 9.5 km s^^ RV difference at 



{l,b) = (0°, —10°). The result changed to 5.2 ± 9.5 km s^^ when the authors accounted 
for probabilities that each given star belongs to either RC. The small change of the result 
is surprising to us (see Sec. [5]). Neither of the above mentioned studies searched for the 
changes in relative kinematics of the two arms on the scales smaller than 5°. 
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There is one problem inseparably involved in the kinematic studies of the Galactic 
bulge that use the RC stars. Not all the stars of the RC region of the CMD belong to the 
RC because there are underlying red giants th at may be at diffe r ent di stances and thus 



Mao fc Paczyhskil (120021 ) these stars dilute 



posses different mean kinematics. As noted by 
the measured kinematic differences between the brighter and fainter RC stars. There are 
two ways one can cope with this fact when comparing the observations with predictions 
based on Galactic models. The first is to report the raw measured properties and compare 
then a with model pre dictions, which are diluted in the same manner as observations (see 



e.g 



NessetaL 



2OI2I ). The disadvantage of this method is that the above mentioned results 
cannot be compared in detail to the models as long as one does not specify from which 
parts of the CMD stars were used in the calculations. The other approach is presented in 
this paper. The probability that a given star belongs to the brighter or fainter arm of the 
X-shaped structure is assigned to each star (based on its extinction-corrected brightness) 
and using this we find the intrinsic properties of the two arms. The result can be directly 
compared to the model predictions, which are based upon the proper motions of stars 
selected based on their distances. It also allows the analytical derivation of the Galactic 
properties. 

We note one more proper motion study in bulge fi eld, which does n ot sho w a double RC 



Clarkson et al 



( I2OO8I ) found proper 



structure. Based on the Hubble Space Telescope data 
motions in a field (/, 6) = (1.25°, —2.65°). These authors assigned photometric distances 
to every star and calculated mean proper motions in distance bins. They concluded that 
the mean circular speed of the bulge stars follows solid body rotation to a cutoff radius 



of 0.3 — 0.4 kpc with a maximum velocit y of 25 km s ^. 



rotation velocity of 62 — 83 km s kpc . 



Clarkson et al. 



his c orresponds to an angular 



(120081 ) also estimated that biases 



in the stellar properties decrease the measured rotation velocity by a factor of two. Thus 
the angular rotation velocity from that study is between 125 and 167 km s~^kpc~^. 
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The aim of this paper is to study the proper motions of the stars belonging to the 
two arms of the X-shaped structure. We use data collected by the OGLE-III survey to 
calculate the proper motions. The third phase of the OGLE project used a camera with 
better resolution than the second phase. In addition the observed sky-area was significantly 
larger and covered some fields with the double RC The results presented here could not 
be achieved using the OGLE-II data. We describe the observations used in this study in 
Section [2l The two following sections present the luminosity function construction and the 
calculation of proper motions of individual stars. Both of these are used in Sec. [5] to find 
the intrinsic proper motion differences of stars in the two arms of the X-shaped structure. 
We end by discussing the implications of our results. 

2. Observations 



The OGLE-III project (IUdalskill2003l ) was conducted with the 1.3-m Warsaw Telescope 
located at the Las Campanas Observatory, Chile. The telescope was equipped with an eight 
CCD chip mosaic camera. The total camera dimension was 8 k x 8 k pixels, with a field 
of view 35' x 35'. The 15/im pixels gave a 0.26" pixel scale, allowing full benefit from the 
excellent seeing conditions. The observations were performed in V and I filters with the 
majority of observations taken in the / band. The V band observations are used here only 
for color information. The exposure time was 120 s. 

The main goal of the OGLE-III project were studies of the microlensing events 
which are effectively found only in the Galactic bulge. The observed bulge area was 
92 square degrees within the range —12° < / < 13° and —7° < 6 < 6°. During the survey 
observations, the sky area monitored decreased with time in order to increase the cadence 
in the fields showing the highest number of microlensing events. The number of epochs 
per field ranged from a dozen or so to ~ 2400, and the time span of observations ranged 
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from about one y ear to eight years . The 



were presented by 



Szymahski et al. 



photometric maps of the OGLE-III bulge fields 



(1201 ll ). We select the fields where the double RC is well 
pronounced and observing coverage was good enough to obtain the precision of the proper 
motions below 1 mas yr~^ in each direction. We were left with three fields named BLG134, 
BLG167 and BLG176, the latter two of which are neighboring. The sky coordinates of the 
fields centers, the number of epochs collected and observing time span are summarized in 
Table [H The sky area corresponding to the one CCD chip is called a subfield. Different 
subfields of a given field are distinguished by the number from 1 to 8 separated by a decimal 
point. 

The astrometric measurements analyzed here are independent of the standard OGLE- 
III reduction, which used Difference Image Analysis. In order to measure the centro i ds of 



the stars, we implernented the effective PSF approach presented by 



and 



Anderson et al 



Anderson fc Kind (120001 1 



(120061 ). The effective PSF is a convolution of the instrumental profile 
with the pixel sensitivity map. Both of the latter are unknown, but they are not needed as 
the stellar centroids can be found using an accurate effective PSF. In our implementation 
the grid on which the effective PSF is defined supersamples the pixel by a factor of four in 
each direction. The spatial variation of the profile is treated by calculating it on a 2 x 4 



Table 1: Characteristics of the observed fields. 



field name 


R.A. 


Dec. 


/ 


b 


-' ' epoch 


At 




J2000.0 


J2000.0 


'o' 


o" 




yrs 


BLG134 


17^57!"6 


-34°12' 


-3.24 


-4.88 


326 


4.6 


BLG167 


IS^OSf'S 


-srso' 


-0.56 


-4.80 


360 


4.6 


BLG176 


18^06?^4 


-31°15' 


0.23 


-5.00 


355 


4.5 



grid for each CCD chip individually. The effective PSF for a given place of the CCD chip is 
evaluated by a linear interpolation of the four nearest grid profiles. 

3. Luminosity function 

In order to assign each given star a probability that it belongs to the brighter or 
fainter arm we der edden each star individually using the interpolated extinction maps of 



Nataf et al. 



( 120121 ). The E{V — I) reddening varies between 0.53 mag and 0.81 mag for 
BLG134 and between 0.62 mag and 1.29 mag for the other two fields. Figure [1] presents a 
sample CMD after the extinction correction was applied (Jq and {V — /)o are extinction 
corrected values of / band brightness and {V — I) color). We use a two-dimensional 
extinction map, so the foreground disk stars are artificially shifted toward bluer colors and 
brighter magnitudes. 

The luminosity function of the bulge red giants is constructed by taking into account 
the stars that are redder than the gray line in Fig. [H This line changes its slope at 
Iq = 14 mag to exclude the disk RC stars, which are abundant in the CMD region 
{V — /)o ~ 0.8 mag and Iq < 13.5 mag. We argue that this selection of stars removes the 
great majority of the foreground disk red giants. The nearby red giants are shifted toward 
bluer colors by extinction correction and thus ar e not included in o ur selection box. We 



20031 ). Additional cleaning 



verified this using the Besangon Galactic model (iRobin et al.l 
of the sample is conducted later on, when the objects with the highest proper motions are 
removed. A separate luminosity function was constructed for each field. All are shown in 
three panels of Fig. [21 We found that calculation of a separate luminosity function of each 
OGLE field is close to the optimum between obtaining good statistics and characterizing 
the spatial changes in the luminosity function. 
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Fig. 1. — Example of a derredened CMD (for the subfield BLG175.6). Stars lying to the 
right of the gray dashed line are used in construction of the luminosity function. Two longer 
horizontal lines mark the brightnesses of the two RCs. The two shorter horizontal lines mark 
the brightnesses of the two RGBBs. 
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The luminosity function for each field is fitted with an analytical functio n that takes 
into account both arms of the X-shaped structure. In previous studies (e.g. 



Natafetal 



20111 ) simpler functions were used. The number of stars in each arm is a sum of three 



components: red giant branch (represented by an exponential fun ction), RC (r e presented by 



a Ga ussian) and the red giant branch bump (RGBB henceforth; 



Gallart 



1998 



Nataf et al 



2011 



; represented by a Gaussian). We take into account the RGBB and ignore the 



asymptotic giant branch bump which was also found in the Galactic bulge (JNataf et al. 



20111 ) because of two major differences between the two structures. First, the number 
counts of the asymptotic giant branch bump are at least an order of magnitude smaller. 
Second, the absolute brightness difference between the RGBB and the RC is significantly 
smaller than for the asymptotic giant branch bump (0.737 mag compared to 1.06 mag). 
The function that we fit has sixteen parameters (two for each of the exponentials and 
three for each Gaussian), but we make several assumptions to reduce the number of fitted 
parameters. The slopes of the exponential components were kept the same. The number 
of the red giants at brightness of the RC was also kept the same. The dispersions of 
the Gaussians representing the RGBB were the same as the dispersions of the Gaussian 
representing RC in each of the arms. The brightness difference between RGB B and RC as 



well a s the number ratio of RGBB to RC stars were fixed at values found by 



Nataf et al. 



( 120131 ) of 0.737 mag and 0.201, respectively. These parameters depend on the metallicity 
and thus their true values in the analyzed fields may be different than assumed. We made 
two exceptions to the fitting procedure described above, which significantly improved our 
fits. First, in the field BLG134 the contribution of the fainter RGBB was neglected. Second, 
in the field BLG176 the number of the RGBB stars relative to the RC stars in the fainter 
arm was a free parameter with the best fitting value of 0.075 ± 0.063, compared to 0.201 in 
the standard fit. The fits had eight or nine free parameters and between 81 and 101 data 
points in the luminosity function. We adjust the interval of Jq in which the fit is performed 
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for each field separately. The difference in the extinction-corrected brightness of the two 
RCs are 0.353 ±0.034 mag, 0.365 ±0.051 mag, and 0.387 ±0.046 mag for BLG134, BLG167, 
and BLG176, respectively. 

Each star was assigned a probability that it belongs to the brighter (ph^Io)) or the 
fainter (pj(/o)) arm of the X-shaped structure. This probability is equal to the ratio of 
number of stars from a given arm to the total number of stars, both of which are the values 
of the fitted function at the extinction-corrected brightness of the star. The plots of Pb{h) 
are shown on the three panels of Fig. [21 

4. Calculation of proper motions 

The crucial point in the calculation of high accuracy proper motions is finding an 
adequate transformation of grids. We note that the OGLE-III observing strategy was not 
optimized for precise astrometry: there was typically only one observation of a given field 
taken per night, the seeing FWHM varied up to a factor of 2.5, exposures were taken at 
different airmasses, and the airmass could change significantly between any two consecutive 
exposures. These issues prevent us from using a detailed ge ometric correcti on common 



Bellini fc BedinI (e.g. 



2010|) 



for all of the subfields of the given exposure as was done by 

Instead, we derived the transformation grids for each subfield individually. This enforced 

the proper motion zero points to be different for each subfield. 

The cross-matching of the star c atalogs from each im age was done using a simplified 



version of the algorithm presented by 



Pal fc Bakod (120061 ). Second order polynomial grids 



were calculated in order to transform the measured positions to a common frame. The 
stellar centroids from all the frames of a given subfield were transformed to this common 
frame. After that, we fit a model for each star that took into account the coordinates for 
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Fig. 2. — Luminosity functions and derived probabilities as a function of extinction-corrected 
Jo brightness. The odd panels (counting from top) present the luminosity functions (gray 
points) and corresponding fits (solid line). The dotted and dashed lines present the contri- 
bution of the brighter and fainter arm, respectively. The even panels show the probability 
that the star of given brightness belongs to the brighter (dotted line) or fainter(dashed line) 
arm of the X-shaped structure. 
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a selected epoch, the proper motion and the differential refraction coefficient. The value 
of differential refraction shift is calculated by multiplying the coefficient by the tangent of 
the zenith distance. The zero point of the proper motions is set by calculating the mean 
of proper motions for red giants brighter than Iq = 13.8 mag. Stars with proper motions 
significantly higher than those of bulge stars are not considered. From the proper motion 
of every star the mean proper motion of bulge stars is subtracted. These proper motions 
and differential refraction coefficients are used to calculate the residuals of the stellar 
positions. Those are transformed back from the common frame to the individual frames 
and the measured positions are subsequently corrected for these residuals. Then, new grid 
transformations are found and using them we calculate the mean position of every reference 
star. If these differ more than 0.15 pix from the initial common frame, we remove the star 
from the list of reference objects. The mean positions are used as a new common frame, 
which in principle is corrected for differential refraction. The grid transformations from the 
individual frames to that frame are found and final fits of positions for a selected epoch, 
proper motions and differential refraction coefficients are performed. 

The uncertai nties of the proper motions are found using the bootstrap method 



(^Press et al. 



19921 ). From a set of exposures of a given field, we draw with replacement a 
subset whose number of elements is equal to the original set. Multiple such subsets are 
drawn and for each of them the procedure of calculating the proper motions described 
above is repeated. The standard deviation of the proper motion of a given star is taken as 
the measurement uncertainty. The proper motion uncertainties are about 0.3 mas yr~^ for 
/ = 14 mag and about 0.5 mas yr^^ for / = 16 mag. It was found that the proper motion 
uncertainties increase near the edges of subfields. This is caused by the less well defined grid 
transformations in these parts of the subfields. We compared the proper motions measured 
for stars present in overlapping parts of the adjacent subfields and found that bootstrap 
estimates are consistent with measured proper motion differences. The raw measurements 
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of proper motions used in this study will be published as a part of the proper motion 
catalog covering the whole OGLE-III bulge area. 

5. Analysis 

The intrinsic (i.e., undiluted by red giants) proper motion difference between the 
brighter and fainter RC can be found using the calculated proper motions and probabilities 
that a given star belongs to either the brighter or fainter arm. Let index i label the 
stars. For each of them we have measured the proper motion in Galactic coordinates: 
/^i,i*, Aij,aJ and corresponding uncertainties ^j^^*, ^j^j,. We know also the extinction-corrected 
brightness Jj^o which gives the probabilities Pb{Ii,o) and pf{Iifi). For the fainter arm the 
average (dispersion) of the proper motions in Galactic longitude is denoted /if,/^ (c^f,;*), 
while for Galactic latitude the corresponding value is /if^b (cf,^)- For the brighter arm the 
corresponding symbols have first index changed from f to b. The likelihood function for a 
single star (£j) is defined as: 






.,bJ 



: exp ' 



2 



2 , f 2 



2-A/Kl^T^!Tfe^Ky V 2KV + e/,J 2ial, + i 



The product of likelihoods for all the stars in a given subfield is the function we 
maximize using the Markov Chain Monte Carlo (MCMC). A separate chain was run for each 
subfield. After trial and error we choose the same interval of Jq brightness between 14 mag 

*The fiu = fiicosb and fib are the proper motions in Galactic coordinate system with 
longitude value corrected for scale changes. 
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and 15 mag. In this range both RCs are prominent, and that is where the most information 
on the undiluted proper motions comes from. For the stars brighter than Jq = 14 mag 
the probabihties are poorly constrained, and the number of brighter stars is smaller. For 
the stars fainter than Iq = 15 mag, the RGBBs significantly contribute to the luminosity 
function and the analytical fits presented in Sec. [3] are slightly poorer. For Iq ~ 15.7 mag, 
the main sequence disk stars start to contribute significantly to the luminosity function. 
The mean value of the effective number of stars in the brighter (fainter) arm i.e., YliPbi^ifi) 
(EiP/(^i,o)) per subfield is 1093 (1499) in the field BLG134. In the field BLG167 the 
corresponding value is 1723 (669), and for BLG176 it is 1416 (771). The values in individual 
subfields do not differ by more than 10% except the subfields BLG167.5 (15% larger than 
the mean) and BLG176.5 (13% larger than the mean). 

In order to illustrate the reliability of our MCMCs, we overlay the Pb(/o) and /Xj,/* 
averaged in 0.1 mag wide /q bins. The panels of Fig. [3] present such plots for four 
sample subfields. The agreement between the measured proper motions and the ph for 
Jo between 14 mag and 15 mag proves the consistency of our procedures of extinction 
correction, luminosity function construction, fitting analytical model to it, and proper 
motion calculation. 

The results of the MCMC are presented in Table O We do not state values of /ib,«+, /^b,6, 
fif^i^,, and yUf.6 as their zero points may change in different subfields. Instead we present their 
differences i.e., A/x^ = /^b,i* — A'-f,/* and A/i^ = /it,,;, — Aif,b, which carry all the astrophysical 
information. 

We note that the differences of /i;^ at the brightness of the two RCs are up to 
1.2 mas yr~^ in the subfields presented in Fig. [31 The values of A/i;^, which are corrected 
for dilution by red giants, for these subfields are about twice larger. Such a large ratio is 
caused by the fact that even for the brightness of the brighter RC some contribution of the 
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14 15 

/„ [mog] 



14 15 



Fig. 3. — Comparison of /i;* (binned in overlapping 0.2 mag wide bins which are 0.1 mag 
apart) and probabihties derived as a function of the dereddened brightness. Each panel 
presents a different subfield. The dashed vertical lines show the range 14 mag < Jq < 15 mag 
that we use in the fit. 
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faine r arm is seen and v ice ve rsa. Among the previous comparisons of two arms kinematics, 
only 



Uttenthaler et al. 



( 120121 ) accounted for the probabihties that measured stars belong to 



the either RC In contrast to us they found a very small change of the kinematic properties. 
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AU analyzed fields are close to 6 = —5°. The range of Galactic longitudes is significantly 
larger. The fields BLG167 and BLG176 span from I = -0.9° to / = 0.6°. The field BLG134 
is located around 2° from these two and is therefore discussed separately. 

Figure m presents the proper motion differences A/i^^ and Afib vs. Galactic longitude. 



We are analyzing 



:he fields close to / = 0°, so the pur ely cylindrical rotation identified 



from RV surveys (JKunder et al. 



2012 



Ness et al. 



20131 ) should result in constant values 



of Afii^ and in Afib being close to zero. However, this expectation actually only holds 
for negative /. For positive /, both A/i/^ and A/i^ are linear functions of /. We fit a five 
parameter phenomenological model to both Afii^{l) and AyU^lO- The subfield BLG176.2 is 
not included in this fit since A/x;^ and Afib derived in this field are the most outlying points 
on both panels of Fig. HI For / < /break we assume 



, const 



AfiUl) = Afil 
Afibil) = A/^r^* 
and for / ^ /break we assume 

Afiuil) = A/ir^* + ^^ ■ (/ - /break) 
Aflbil) = A/ir^* + ^|P ■ (l - ^break) • 

The fit resulted in /break value of —0.114° ± 0.060°, i.e.,2cr away from 0°. The 
parameters for the longitudinal proper motion are A/i™"'** = 2.04 ± 0.07 mas yr^^ 
and —jf^ = 1-23 ± 0.32 mas yr^^ deg^^. For the latitudinal proper motion we obtain 
^^const ^ _o.i4±o.08 mas yr^^ and ^ = -3.27 ±0.43 mas yr^^ deg-\ The fit resulted 
in xVdof = 42.1/25. 
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Table 2. Proper motion statistics of two arm of the X-shaped structure 



subficld 



in 



Art. 



0"b,i* 



O'id* 



^IJ-b 



"■b.ii 



o"f,f) 



BLG134.1 
BLG134.2 
BLG134.3 
BLG134.4 
BLG134.5 
BLG134.6 
BLG134.7 
BLG134.8 
BLG167.1 
BLG167.2 
BLG167.3 
BLG167.4 
BLG167.5 
BLG167.6 
BLG167.7 
BLG167.8 
BLG176.1 
BLG176.2 
BLG176.3 
BLG176.4 
BLG176.5 
BLG176.6 
BLG176.7 
BLG176.8 



-3.358 

-3.226 

-3.095 

-2.963 

-3.108 

-3.242 

-3.375 

-3.507 

-0.680 

-0.547 

-0.415 

-0.281 

-0.428 

-0.560 

-0.696 

-0.830 

0.104 

0.240 

0.373 

0.507 

0.366 

0.230 

0.097 

-0.037 



-5.129 
-5.053 
-4.976 
-4.900 



1.95 ±0.23 2.90 ±0.10 1.91 ± 0.07 
1.71 ±0.23 3.04 ±0.10 2.21 ± 0.08 



1.72 ±0.25 2.97 ±0.11 1.99 ± 0.08 

1.81 ±0.19 2.87 ±0.09 1.88 ± 0.07 

-4.646 3.25 ±0.20 2.84 ±0.13 2.28 ± 0.07 

-4.720 2.25 ±0.24 2.82 ±0.12 1.84 ± 0.07 

-4.795 1.89 ±0.27 3.09 ±0.10 2.14 ±0.07 

-4.872 3.37 ±0.18 2.38 ±0.10 2.31 ± 0.08 



-5.039 
-4.967 
-4.898 
-4.821 
-4.556 
-4.633 
-4.706 



2.05 ±0.18 2.73 ±0.07 1.45 ±0.11 

2.40 ±0.21 2.73 ±0.07 1.83 ±0.11 

1.89 ±0.19 2.79 ±0.07 1.46 ±0.10 

2.16 ±0.21 2.64 ±0.06 1.58 ± 0.09 

2.29 ±0.20 2.82 ±0.07 1.78 ±0.11 

2.24 ±0.20 2.80 ±0.07 1.72 ± 0.09 

1.88 ±0.25 2.87 ±0.08 1.76 ±0.14 



-4.781 2.00 ±0.20 2.49 ± 0.07 1.75 ±0.13 



-5.245 1.97 ±0.22 2.84 ± 0.08 



1.80 ±0.10 

3.33 ±0.10 1.97 ±0.10 

2.69 ±0.10 1.90 ±0.09 

2.94 ±0.19 2.62 ±0.08 2.02 ± 0.09 

2.84 ±0.23 2.77 ±0.09 2.21 ±0.10 

2.31 ±0.17 2.90 ±0.07 1.84 ± 0.09 

2.14 ±0.27 2.82 ±0.09 2.07 ±0.11 

1.79 ±0.23 2.85 ±0.08 1.75 ±0.10 



-5.176 3.42 ±0.21 

-5.104 2.62 ±0.22 

-5.029 

-4.772 

-4.847 

-4.921 

-4.994 



-0.11 ±0.19 2.68 ±0.09 1.85 ± 0.08 

-0.08 ±0.21 2.97 ±0.10 1.87 ±0.08 

0.23 ±0.20 2.67 ±0.11 2.01 ± 0.08 

-0.18 ±0.18 2.78 ±0.08 1.78 ± 0.06 

-1.93 ±0.17 2.84 ±0.11 2.27 ±0.07 

-0.14 ±0.21 2.78 ±0.11 2.13 ±0.09 

-0.65 ±0.17 2.85 ±0.08 2.00 ± 0.07 

-2.74 ±0.19 2.26 ±0.10 2.62 ± 0.09 

0.40 ±0.18 2.63 ±0.07 1.76 ±0.12 

-0.07 ±0.25 2.58 ±0.07 2.02 ±0.13 

-0.34 ±0.20 2.54 ±0.06 1.83 ±0.13 

-0.04 ±0.20 2.64 ±0.07 1.91 ±0.13 

-0.38 ±0.22 2.76 ±0.07 2.01 ±0.14 

0.06 ±0.23 2.87 ±0.06 2.03 ±0.14 

-0.24 ±0.26 2.63 ±0.06 2.48 ±0.13 

-0.66 ±0.28 3.22 ±0.08 2.49 ±0.17 

-0.60 ±0.24 2.79 ±0.07 1.94 ±0.12 

-3.09 ±0.24 3.20 ±0.10 2.23 ±0.12 

-1.52 ±0.21 2.61 ±0.08 1.84 ±0.10 

-1.97 ±0.20 2.61 ±0.08 2.30 ±0.10 

-1.76 ±0.20 2.69 ±0.08 1.96 ±0.11 

-1.66 ±0.21 2.61 ±0.07 2.15 ±0.10 

-1.20 ±0.21 2.95 ±0.07 1.93 ±0.11 

-0.24 ±0.26 3.28 ±0.10 2.30 ±0.14 



Note. — First three columns give the name of the subfield and the Galactic coordinates of its center. All proper motion 
differences and disperions are in mas yr^^. 
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[deg] 



Fig. 4. — Proper motion differences of two arms of tlie X-sliaped structure as a function of 
tlie Galactic longitude. The upper panel shows the longitudinal proper motion difference 
while the lower one presents the latitudinal proper motion difference. The lines are fits with 
common longitude of the break point in both panels. 
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6. Discussion 



The most striking result of Fig. H] is the presence of significant changes of the proper 
motion differences with the Galactic longitude for / > —0.1°. We interpret this as a 
signature of asymmetric streaming motions of stars in the X-shaped structure. The fact 
that the mean motion of the stars belonging to a certain structure is different from the 
mean motion of this structure is known e.g., in the case of Galactic spiral arms. The 
measurements presented in this paper are insufficient to determine in which of the arms the 
streaming motion occurs. None of the published models of the X-shaped structure gives 
predictions of the streaming motions that we report. 

The value of A/i^^"''* can be used to constrain the angular velocity of the Galactic bar. 
For small values of distance modulus differences of two RCs A/rc their relative distances 
can be approximated by Ji12__rc_ ^jj^g ]^g^j. angular velocity depends on the proper motion 
difference measured and the longitudinal proper motion of Sgr A*: 

5^^const 



n 



Z^SgrA*,;* + 



InlOA/^c 



Substituting the A/x™'''^*, AIrc i n BLG167 of -0.3 6 5 ± .051 mag and /iggr a*m 



Reid fc Brunthalerl toOJ ) we obtain Q = -18.5 ± 



—6.38 mas yr ^ as measured by 

1. 8 mas yr~^ = — 8 7.9 ± 8.2 km s^^ kpc^lJ. This result is marginally consistent with results 



of 



Clarkson et al 



( 120081 ) without correcting for biases in their analysis. If we assume that 
the Galactic angular velocity is constant up to some Galactic radius R' and for larger radii 
the circular spee d is the same as the local circular velocity of the Sun Vc = 238 ± 9 km s~^ 



( jSchonrich 



2OI2I ) then we can estimate R' = Vc/ \ fl\ = 2.71 ± 0.28 kpc. The bar angle 



relative to Sun-Galactic center line is ~ 30° (e.g., 



Cao et al 



20131). Thus we can estimate 



that the bar points that are 2.71 ± 0.28 kpc away from the Galactic center are observed at 



^Negative sign indicates the same direction of bulge rotation as the Sun's orbital motion. 
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l = 12.87 ± 2.4° and / = —7.27 ± 0.44° (these uncertainties do not take into account the 
uncertainty of the bar angle). 

The proper motion differences in the field BLG134 are similar to those found in field 
BLG167. The two most outlying subfields are BLG134.8 and BLG134.5. The former has 
the smallest longitude among the subfields analyzed here, while the latter is the closest to 
the Galactic disk among the BLG134 subfields. Thus, these discrepancies are most likely 
caused by significantly smaller contribution of the brighter arm in these subfields. 

The ratio of proper motion dispersions in both arms along the latitudinal direction 
can be used as an independent estimate of the distance ratio of the two arms if the same 
tangential velocity distribution is assumed in each arm. We note that this may not be the 
case since in a given sight line the two arms are at the different distances from the Galactic 
disk. 

T here are several studies devoted to the proper motion dispersions in the Galactic 



bulge (jSpaenhauer et al. 



1992; 



Kuijken fc Rich 



2002; 



Vieira et al. 



20071 and references 



therein). These dispersions were found to be anisotropic with cru/cxb ~ 1-2. We do not 
confirm this finding when the two arms are treated separately. Analysis of the posterior 
probability distribution of the MCMC revealed that only in five cases did the ratios differ 
significantly from unity. In each of these cases the longitudinal proper motion dipsersion is 
smaller than the latitudinal one. The value of the dispersion ratio in the brighter arm i.e. 
(ThM/(^h,b is 0.868 ± 0.032 in BLG176.8 and 0.774 ± 0.025 in BLG167.8, while for the fainter 
arm at^ijat^b is 0.763 ± 0.071, 0.707 ± 0.072, and 0.709 ± 0.063 for BLG176.8, BLG167.8, 
and BLG167.7, respectively. 
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